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ABSTRACT 

We calculate the prompt x-ray emission as a function of viewing angle for beamed 
Gamma-Ray Burst (GRB) sources. Prompt x-rays are inevitable due to the less 
highly blueshifted photons emitted at angles ;> I/7 relative to the beam symmetry 
axis, where 7 is the expansion Lorentz factor. The observed flux depends on the 
combinations "fAO and r y0 v , where AO is the beaming angle and V is the viewing angle. 
We use the observed source counts of 7-ray-selected GRBs to predict the minimum 
detection rate of prompt x-ray bursts as a function of limiting sensitivity. We compare 
■ our predictions with the results from the Ariel V catalog of fast x-ray transients, and 

. find that Ariel's sensitivity is not great enough to place significant constraints on 7 

and AO. We estimate that a detector with fluence limit ~ 10~ 7 erg cm -2 in the 2-10 
I keV channel will be necessary to distinguish between geometries. Because the x-ray 

emission is simultaneous with the GRB emission, our predicted constraints do not 
involve any model assumptions about the emission physics but simply follow from 
CIh 1 special-relativistic considerations. 

6 

^ ' Subject headings: gamma rays: bursts 
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1. Introduction 

With the 1996 launch of the BeppoSAX satellite (Boella et al. 1997), it became possible 
to localize Gamma-Ray Burst (GRB) sources to arcminute accuracy within a few hours after 
their detection. Such localizations were followed quickly by the discovery of delayed counterparts 
("afterglows") in the x-ray (Costa et al. 1997), optical (van Paradijs et al. 1997), and radio 
(Frail et al. 1997) regimes of the spectrum. Subsequent identifications of spectral features in the 
afterglows and host galaxies of four GRBs led to the determination of source redshifts z = 0.835 
for GRB 970508 (Bloom et al. 1998), z = 0.966 for GRB 980703 (Djorgovski et al. 1998), 
z = 3.418 for GRB 971214 (Kulkarni et al. 1998), and z = 1.6 for GRB 990123 (Bloom et 
al. 1999). Thus the distance scale was firmly established as cosmological, and the energy scale 
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deduced from the observed fluences was very large, ^ 10 51 -10 54 ergs for isotropically-emitting 
sources. The standard picture which naturally follows from such an enormous energy release is 
a relativistically-expanding fireball with a Lorentz-factor 7 ;> 10 2 , which produces the observed 
gamma-rays with the formation of either internal (Paczyhski & Xu 1994; Rees & Meszaros 1994; 
Sari & Piran 1997; Pilla & Loeb 1998) or external (e.g. Meszaros & Rees 1993) shocks. Examples 
where the inferred energy release for isotropic emission exceeds the expected radiative energy 
supply from compact stellar-mass objects (Kulkarni et al. 1998; Fruchter 1998; Kulkarni et al. 
1999), suggest the existence of beaming. 

The degree of beaming in GRB sources is still poorly constrained. The expansion could be 
taking place in a jet whose opening half-angle is anywhere from C I/7 (strongly beamed) to tt 
(isotropic expansion). Many of the current models for GRB progenitors involve systems which 
have a rotation axis (e.g., Paczyhski 1998; Fryer et al. 1998; Popham et al. 1998; Woosley 1998; 
Meszaros et al. 1998), possibly leading to a preferred expansion direction for the explosion. With 
beaming, the total energy necessary to produce the observed fluxes will be reduced by a factor 
AS7/47T, where A$7 ^ I/7 2 is the solid angle subtended by the GRB radiation. In addition, 
the total event rate must be greater by a factor 4n/ACl. As the Lorentz factor of the emitting 
material declines with time, it is likely that the delayed emission at longer wavelengths will be less 
relativistically beamed than the jet collimation angle, and that detection rates of afterglows (after 
correcting for different detector sensitivities) will be correspondingly larger (Rhoads 1997; Perna 
& Loeb 1998). Grindlay (1999) has investigated this possibility by comparing the discovery rate of 
GRBs by the BATSE experiment with that of fast x-ray transient sources (those which exceeded 
the flux threshold for only ~ 1 orbit) by the Ariel V satellite (Pye &: McHardy 1983) two decades 
ago; he found that the rates were roughly consistent with no differential beaming between the 
7-ray and the x-ray emission. The naive interpretation of this constraint is, however, sensitive to 
the assumption that each GRB is followed by a bright X-ray afterglow. In this paper, we try to 
avoid this assumption and consider instead the inevitable (simultaneous) emission of X-rays at an 
angle relative to the GRB beam, purely due to special relativistic considerations. 

If the 7-ray emission is beamed, it is unavoidable that there will be some photons which 
are detectable even when the observer is not within the 7-ray emitting beam. These are the 
photons emitted at large angles to the expansion direction by material near the edge of the jet. 
Such photons will not generally be detectable as 7-rays; since they are less highly blueshifted 
than forward-emitted photons, they will be seen in the x-ray band. Thus, there should be some 
minimum degree of differential beaming between the 7-rays and the prompt x-rays. The degree 
to which this will affect the x-ray detection rate depends on the expansion Lorentz factor 7 and 
the beaming half-angle A9; comparison of the predicted rate with that observed by Ariel V can 
yield constraints on 7 and A6>. 

In this paper, we use the BATSE source counts and a simple geometrical model for the GRB 
sources to predict the prompt x-ray detection rate as a function of limiting flux. In §2, we derive 
the expression for the flux as a function of viewing angle for a given geometry. In §3, we calculate 
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the x-ray detection rate for various values of 7 and AO, and compare the results with the rate 
observed by Ariel V. In §4, we present our conclusions. 



2. Calculation of the X Ray Flux 

2.1. Emission from Optically Thin, Relativistically Moving Material: General Case 

It will be instructive first to calculate the observed flux for the general case of emission 
from optically thin, relativistically-moving material. In doing so, we closely follow the analysis 
of Granot, Piran, & Sari (1998). Suppose we have some optically-thin material with rest- 
frame emissivity j' measured in erg s _1 cm~ 3 Hz -1 sr _1 . Following from considerations of 
Lorentz-invariance (Rybicki & Lightman 1979), the lab-frame emissivity per frequency v is given 
by 

j = hd. m 

3v 7 2 U-/^) 2 ' {) 

where 7 = (1 — /3 2 ) -1 / 2 is the local Lorentz factor, /j, v is the cosine of the angle between the local 
velocity and the direction to the detector, and v' = i"f(l — PfJ, v ). In general, j' u , is a function of 
position, time, frequency, and direction; we are only concerned with photons which are emitted in 
the direction £l' d toward the detector. 

Let us use a spherical coordinate system r = (r, #,</>), where the coordinates are measured 
in the lab frame; let the 9 = axis (z-axis) point toward the detector. Suppose the detector 
is at a distance D from the source at the origin (see Figure 1). Furthermore, let a be the angle 
that a given ray makes with the normal to the detector. The flux at the detector, measured in 
erg s _1 cm~ 2 Hz -1 , is then given by 

I v (a,(f)) d£l = / d(f> da sina cosq I v (a,(f)), (2) 
Jo JO 

where I u (a,(p) is the intensity along a ray incident on the detector in direction (ot,4>). We will 
only consider situations where the source is very far away, so that all rays come in very nearly 
normal to the detector. Note that in this case, a <C 1, so sina w a and cos a ps 1; furthermore, 
the only significant contribution to the integral comes from a ^ L/D « 1, where L is the size of 
the source. In the approximation that a€l, simple geometry gives 



r v j 

a = —sm9 = —^1- fi 2 , (3) 



where we have defined fi = cos 6. The equation of a ray incident at angle a is thus 

aD 



smV 

so the change in r corresponding to an increment in 9 along the ray is given by 



(4) 



dr = -aD^-dO. (5) 
sin^ 9 
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Therefore, the line element along a ray is given by 

ds = {r^ + dr^ = J ^ W2 d,. (6) 

For an optically-thin medium, the contribution to the intensity from a ray segment ds is 
dl u = j u ds. Note that if this dl u is received at the detector at time T, then it relates to the 
emissivity j v at an earlier time t, due to the light -travel time. Let T = be the arrival time at the 
detector of a photon emitted at the origin at t = 0; inspection of the geometry in Figure 1 gives 

t = T+ r -^. (7) 

c 

Combining equations ([!]), @, (a), and (^) gives 

Here, a m = L/D where L is the projected size of the emitting region; £l' d is the direction toward 
the detector, as measured in the rest frame; and r = (r,9,<fi) = {aD/yl — /i 2 , cos -1 /x, <j>) is the 
lab-frame position vector. 

Equation (|8|) is quite general; if we can specify the rest-frame emissivity (tt' d , r, t), perhaps 
from some physical model of the radiative processes taking place, then we can calculate the flux at 
the detector. Note that in general, 7 (and therefore (3) will be functions of r and t; this dependence 
can often be derived from hydrodynamical considerations. 

In the case of radial expansion, fj, v = jx, and we may switch integration variables from \i to u', 
and rewrite equation (||) as 

vD f 2* f a m 9 rnd+Vdj/ faia^T +(!-£) fc] 



F V (T) = -^r \ dcf> a 2 da / ^ ^±£21 . (9) 

K> iPJo V Jo U' 2 {1 - [1(1 - ^)]2}3/2 ^ 



2.2. GRB Emission from a Conical Section of an Expanding Spherical Shell 

Let us assume that the expanding material is confined to a cone of opening half-angle A#, 
and that the emission takes place in a thin shell which moves radially outward at the expansion 
speed /3c, with associated Lorentz factor 7. Realistically, based on hydrodynamical considerations, 
the shell should be of radial thickness ~ 0.1/?ci, but our results are not changed significantly by 
assuming that it is infinitesimally thin. 

The observed spectral shape of GRBs is well-reproduced by a broken power law. Usually, the 
two segments with different power-law indices are joined smoothly with an exponential transition 
(Band et al. 1993). Recently, Ryde (1998) devised a single analytical function which also produces 
a smoothly-broken power law. We adopt this functional form for the frequency dependence of the 
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rest-frame emissivity, since it yields a spectral shape similar to that observed (with the same form 
but with a blue-shifted break frequency and a broader power-law transition). We assume typical 
power-law indices, F' v oc (z/)° for low frequencies and F' v oc (z/) _1 for high frequencies (Band et 
al. 1993). In the case where the axis of the jet points directly at the observer, we have 



j' v ,(r, t) = A(t) 5(r - (3d) H(A8 - 8) 



(10) 



Here, A(t) is the time-dependent normalization, and H(x) is the Heaviside step function, v' is the 
transition frequency, and r) = ln(l + Az/ /vq), where Az/ is the width of the power-law transition. 
In order to achieve agreement with the results of Strohmayer et al. (1998) for the average ratio 
of x-ray to 7-ray flux of Gmga-detected bursts, F(2-10 keV)/F(50-300 keV) = 0.24, we require 
that the observed spectral break occur at his® = 18 keV. This is low compared to the typical break 
energy ~ 150 keV for BATSE bursts (Band et al. 1993), since Ginga was sensitive to lower-energy 
photons. Note that even the BATSE bursts span a range of break energies (10-10 3 keV); and since 
Ginga is one of the few satellites to have probed the prompt x-ray emission from GRBs, we adopt 
the Ginga value. Given that for a typical GRB expanding with Lorentz factor 7^1 toward the 
observer, the rest-frame frequency corresponding to vq is v' Q ps 1^/(27), we use hi/ = 9 1 keV. 
The width Az/ is unimportant for our analysis, since it does not appreciably affect the fluxes we 
calculate; we assume hAv' = 1 keV. The observed GRB spectral shape predicted from equation 
( |Tfj| ) will have the same form, but with the break frequency z^o blueshifted, and the transition 
width Av broadened. Of course, in general, v'§ will be a function of time; for simplicity, we assume 
that it is a constant (i.e. constant spectral shape at the source). Note also that z/jV oc const for 
z/ v'q, which means that we need to introduce a high-end cutoff frequency u' c so that j' v i will not 
possess a logarithmic ultraviolet divergence. We set this upper cutoff to correspond to an observed 
photon energy of 5 MeV. 



Recalling that v' = v~f(l — (3 cos 8), we may combine equations (|9|) and ( |i~0| ) to obtain 



2i:vD 



a 2 da 



fin(l-0 cos A0) du i 



,/2 



1 + 



{i 



1 



(ii) 



where a m = L/D is the angular size of the source when A(t) effectively cuts off, and equation (|3|) 
gives the relation between r and a with \x = /3 _1 [1 — (y' jv^)\. Note that equation (11) implies that 
the observed effective duration of the event is inversely proportional to the observed frequency. 
This is due to the longer travel time for the less-blueshifted photons coming in from off the line of 
sight. 



Now, recall the following property of the Dirac delta function: for f{y) monotonic, 

df- 1 



s[f(y)} 



dy 



<% - yo) 



(12) 
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where yo is the zero of f{y). This allows us to transform the delta function in equation ( |Tl|) to a 
delta function in a. Thus, we can perform the a-integration in equation (|TTj) to obtain 



F U {T) = 2npY — 



dry 



D J 



1/7(1-/3 cos AO) 
1/7(1-/9) 



775 ' 



(13) 



The integral may be computed numerically for any observed frequency v and time T. The 
factor (cT/D) 2 reflects the increase in angular size of the shell as it expands. The case AO = tt 
corresponds to spherical expansion. For simplicity, we adopt the form A(t) = Aoexp(— t/r); this 
produces a single-pulse lightcurve, with total fluence 



S v 



F v (T)dT 



A-kPA TU (CT\ 2 ^7(1-/3 cos A0) 



7 



D 



1/7(1-/3) 



1 + 



dv' 

IF' 



(14) 



In the case where the line of sight to the observer (0 = 0) lies entirely outside the cone, the 
geometry gets a bit more complicated. Let V be the angle between the direction to the detector 
and the axis of the emission cone (see Figure 1). Without loss of generality, we may also place the 
axis of the cone at <fi = 0. In this case, the law of cosines for spherical triangles (e.g., Peebles 1993) 
yields the following equation in (0, (j)) coordinates for the boundary of the cone: 



cos ( 



, cos + sin 9 V sin 9 cos 1 



cos A9. 



(15) 



This translates into a more complicated angular dependence in equation ((T^). We now have that 
j' v , is non-zero for \6 — 9 V \ < A9; for a given 9, the limits on (j) are determined by solving equation 
(|i"5|). Thus, on the right-hand side of equation(|ic|), we have instead of H(A9 — 9): 

cos A# — cos 9 V cos I 



H(A9- 
We then obtain for the fluence 



9J) H 



COS ( 



sin Vv sin ( 



S v 



4pA ru /cr\ 2 

1 \d) 



vy[l-/3cos(6 v +A6)] 
i/7[l-j8cos(0«-A0)] 



A<t>{v') 



1 + 



dv' 



(16) 



(17) 



where solving from equation (15) with cos# = n = (3 1 (1 — u'/'yu) yields 

' cos A0 - f3~ l {l - v' hv) cos 9 V 1 



cos 



sin 9 V y/1 - - v'hvW 



(18) 



This is the same result as equation fll4|), except that the limits of z/'-integration have changed, 
and the (^-integration now gives 2Ac/> instead of 2tt. 

In the intermediate case, where the observer is inside the cone but off-axis, the flux is given 
by the sum of the right-hand sides of equations (|l^) and (|i~7|), with the upper limit of integration 
in (pi) replaced by v^\Y - f3cos(9 v - AO)]. 
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3. Calculation of Prompt X-Ray Source Counts 

3.1. Derivation the X Ray Rate from the Gamma-Ray Rate 

We would like to calculate the expected detection rate of x-ray transients as a function of 
limiting fluence, due to prompt emission from GRBs. In particular, we consider the energy channel 
of 2-10 keV probed by Ariel V (Pye & McHardy 1983), and compute the rate of X-ray transients 
based on the GRB source counts found by BATSE in the 30-2000 keV band (Bloom, Fenimore, & 
in 't Zand 1996). We use fluence instead of flux because the Ariel Fdata have very coarse time 
resolution; the threshold is based on the number of counts detected by the satellite over one orbit 
(~ 100 min = 6 x 10 3 s), and is reported as S x = 2.4 x 10~ 6 erg cm -2 in the 2-10 keV range (Pye 
& McHardy 1983; Grindlay 1999). Our analysis yields the result that the observed duration at 
x-ray frequencies goes roughly like r x ~ r 7 (^ 7 /^ x ), where t 7 is the duration as observed in gamma 
rays, and z^ x and za, are the effective observed frequencies in x-rays and gamma rays respectively. 
Since typically r 7 <J 10 s, we have that r x <; 10 3 s. Thus, the integration time for Ariel Fwas 
sufficiently long for it to detect the entire fluence of the prompt emission. 

If the gamma-ray emission is beamed into a cone of half-angle A8, then we will only detect 
a GRB source in gamma-rays if our line of sight falls within an angle ~ I/7 of the edge of the 
beam. Since the viewing angles 8 V should be distributed isotropically; it follows that the fraction 
of sources whose viewing angles are in the range (81,82) is (cos#i — cos6>2)/2. Thus, the actual 
event rate is enhanced over the detected event rate by a factor 

47r _ 2 

AH ~ 1 -cos[max(A0,l/7)]' 



Note that equations ( |T3[ ) and (o) allow one to calculate the ratio of the Ariel V x-ray fluence 
S x to the BATSE gamma-ray fluence S 7 : 



s x _ ; 2 10 s v {o v )dv 



(20) 



where the frequency v is in units of keV. 

The number of Ariel V x-ray sources with fluences brighter than S x is 

N > SJ " l-co S [J<AM/ 7 )] I?' ~ n(S ^ dS ->- < 21 > 

where n(S^)dSj is the number of sources with fluences in the seen by BATSE, 

and V is the solution for V in equation (|20|), since S x /S^ is a monotonic function of V (cf. Fig. 
2). 
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3.2. Results 

We apply equations ( p0| ) and (plj), with the form of F U (T) derived in §2, to predict the source 
detection rate in the x-ray (2-10 keV) range probed by Ariel V. We find that the results depend 
only on the combinations j9 v and 7A#, as expected from ultra-relativistic beaming. We calculate 
the source counts for four beaming angles: 7A0 = 0.1, 1, 10, and 100. 

In Figure 2, we plot the x-ray fluence as a function of viewing angle 9 V . Clearly, for a given 
event the overall normalization depends on a measurement of, say, the BATSE gamma-ray fluence, 
so we plot the ratio S X (9 V )/S^(0) [cf. Eq. (pO])]. For comparison, we also show the gamma-ray 
fluence as a function of viewing angle, plotted as the ratio 5 7 (^)/5 7 (0). Note that the gamma-ray 
fluence generally falls off to within an angle 7 -1 away from the beam's edge, as expected from 
special relativity. The x-ray fluence decreases somewhat more gradually with viewing angle; this 
effect is significant out to a few times 7 _1 away from the beam's edge. This is due to the relative 
redshifting of photons at larger viewing angles, and is the basis for our assertion that for strongly 
beamed bursts (to within ~ 7 we should see more prompt transient events in x-ray searches 
than in gamma-ray searches. The next step is to apply Equation (^l|) to see how sensitively the 
detection rate should depend on the beaming angle. 

In Figure 3, we plot the source detection rate as a function of limiting flux [cf. Equation 
(|2l|)]. We see that the counts only differ appreciably for limiting fluences fainter than a few times 
10 -7 erg cm -2 . For comparison, we display the Ariel Fdata point, with la Poisson errors. The 
Ariel data are clearly marginally consistent at the la level with any degree of beaming, and 
we would need an experiment at least an order of magnitude more sensitive to be able to place 
meaningful constraints on the beaming. Finally, we note that the results for jA9 ~ 100 agree 
with the log N — log S plots obtained from BATSE (Bloom, Fenimore, &; in 't Zand 1996), once 
one corrects for the frequency band. 

4. Conclusions 

For a given GRB spectral shape, the predicted number of X-ray transients with no GRB 
counterparts depends only on the combination ^yA9, where A9 is the beaming angle and 7 is the 
expansion Lorentz factor of the sources. The excess of X-ray transients in the 2-10 keV band 
becomes significant only at fluences S x <^ 10 -6 erg cm -1 ; therefore, the sample of fast x-ray 
transients detected by the Ariel V satellite in the 1970's (with sensitivity 2.4 x 10~ 6 erg cm -1 ) 
is not deep enough to place model-independent constraints on the beaming angle. It will require 
more sensitive instruments, such as the MAXI all-sky monitor (Kawai et al. 1996) aboard 
the soon-to-be-launched International Space Station, or MOXE on Spectrum-X (in 't Zand, 
Priedhorsky, & Moss 1994), with a comparable or larger sample size, to rule out strong beaming 
at the 2cr-level. In addition, we may hope to use this analysis in the near future to learn more 
about the low-frequency spectral index in GRBs. 
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In contrast with previous limits (e.g., Grindlay 1999; Perna & Loeb 1998), our analysis is not 
sensitive to any model assumptions about the physics of the GRB/afterglow emission, but rely 
only on special relativistic considerations and the observed properties of GRBs. It is inevitable 
that strongly-beamed GRBs will be accompanied by a larger number of detectable prompt x-ray 
events; it is only a matter of pushing the sensitivity limits down by an order of magnitude or so 
before we might expect to observe this excess. 

This work was supported in part by the NASA grants NAG5-7768 and NAG5-7039. 
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Fig. 1. — Geometry for calculating the flux received from a relativistically expanding source. Top 
panel shows the case of a jet aligned with the line of sight; bottom panel shows the case of an 
off- axis jet. 
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Fig. 2. — X-ray fluence S x (solid curves) and gamma-ray fluence 5 7 (dotted curves) as a function 
of viewing angle 9 V , in units of the on-axis 7-ray fluence, 5 7 (0). From left to right, ~/A9 = 0.1, 1, 
10, and 100. 
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Fig. 3. — X-ray (2-10 keV) source detection rate as a function of limiting fluence, for 7A0 = 0.1, 
1, 10, and 100. The data point is the rate of fast x-ray transient detection by Ariel V, which had 
a sensitivity of S x = 2.4 x 10~ 6 erg cm~ 2 in the 2-10 keV channel. The error bar is based on the 
estimated Poisson uncertainty. 



